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Surface structure of quark stars with magnetic fields 

Prashanth Jaikumar* 
Physics Division, Argonne National Laboratory, Argonne IL 60439, USA. 

We investigate the impact of magnetic fields on tfie electron distribution in the elec- 
trosphere of quark stars. For moderately strong magnetic fields B ~ fO^'^G, quantization 
effects are generally weak due to the large number density of electrons at surface, but can 
nevertheless affect the spectral features of quark stars. We outline the main observational 
characteristics of quark stars as determined by their surface emission, and briefly discuss 
their formation in explosive events termed Quark-Novae, which may be connected to the 
r-process. 
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I. INTRODUCTION 



There is renewed interest in the theory and observation of strange quark stars, which are believed to 

contain, or be entirely composed of, deconfined quark matter [1]. An observational confirmation of their 

existence would be conclusive evidence of quark deconfinement at large baryon densities, an expected 

feature of Quantum Chromodynamics (QCD). Furthermore, discovery of a stable bare quark star affirms 

the Bodmer-Terazawa-Witten conjecture [2], viz., that at high enough density, strange quark matter, 

^ ' composed of up, down and strange quarks, is absolutely stable with respect to nuclear matter. This 

0^ . intriguing hypothesis is over three decades old, and bare quark stars are but one possible realization 

C^ ' put forward in the intervening years. Broadly speaking, the field of quark star research follows two 

intertwined paths: (1) efforts to describe the global (mass and radius) and surface structure of such stars 

from microscopies and (2) determination of observational features that distinguish them from their close 

cousins, neutron stars. It is possible that a few of the approximately 1500 neutron stars detected thus 

(*-^ . far are really quark stars [3]. 

r^ ' A positive identification of a quark star would likely require a complement of consistent signals. For 

CIh! example, physical properties such as maximum mass, radius, minimum spin-periods and neutrino/photon 

cooling rates (and hence temperature) all depend on the (unique) equation of state for dense quark 
flj ' matter. As the latter is uncertain, there is a range of allowed values for the above quantities that makes 

(~| , a distinction based on current observational data difficult. We emphasize that a bare strange star is self- 

bound [2,4] while a neutron star requires gravitational forces for its binding. This implies a very different 
equation of state for the two classes of stars; however, in the range of masses (1.35 < M/Mq < 2.4) 
J^ ■ observed to date, the calculated radii are similar {R ^ 10 to 15 km). Perturbative corrections to the 

H ' quark matter equation of state and color superconductivity effects further complicate the issue [5] . ^ 

It was first pointed out by Alcock et al. [11] that the bare quark star surface would be characterized 
by an abrupt termination of the quark density (density gradient ~ 10^®g/cm^) and a less rapid fall- 
off of a charge neutralizing skin of electrons termed the electrosphere (density gradient ~ lO^^g/cm^). 
The latter extends up to 1000 Fermis above the bare quark star surface. The electric field which binds 
these electrons to the surface is large enough to pull electron-positron pairs out of the QED vacuum, 
prompting copious pair-emission which is directed outward by the radiation output by the star [12]. This 
mechanism of emission from the bare surface of a strange quark star, due to both photons and e^e~ 
pairs, can produce luminosities well above the Eddington limit {^ 10^^ erg sec~^) up to a decade after 
the quark star's birth [13]. The spectrum of emitted photons is significantly harder than that from a 
normal cooling neutron star (30 < E/keV < 500 instead of 0.1 < E/keV < 2.5). After the star has 
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We will not elaborate here on the possibility of having superconducting phases of quarks. We refer the interested 
reader to Refs. [6-10] for a consideration of these phases, and the consequences of their existence inside compact 
stars in the present context. 



cooled to temperatures less than T < 6 x lO^K, other microscopic processes of photon emission can 
dominate the luminosity of the star and yield a non-thermal spectrum [14]. This distinctive spectrum 
and temperature evolution, if observed, would constitute an unmistakable signature of a strange quark 
star, but would require a fortuitous detection before the star cools below the observational threshold of 
present-generation space-based laboratories such as INTEGRAL [15]. 

Quark stars may be much more elusive if they have a suspended crust of normal matter up to the 
neutron drip density, that is separated from the quark surface by a microscopic gap but supported by 
the large electric fields at surface [11,16]. Alternatively, quark stars may admit a mixed phase crust 
several meters in thickness and extending up to its surface, if surface and Coulomb costs (determined 
at the microscopic level by the current quark mass of the strange quark) is small enough [17,18]. This 
would obviate the need for large electric fields at surface, smoothen the density gradient at surface, and 
render the previously proffered signals inapplicable. These uncertainties arc not expected to be resolved 
quickly since they involve quantifying the parameters of QCD to better than 5% in a distinctly non- 
perturbative regime. (NB. The density that characterizes quark matter in these objects corresponds to 
a quark chemical potential jiq ^ 500 MeV which does not greatly exceed Aqcd ^ 200 MeV.) 

However, we can certainly aim to quantitatively improve our current understanding of the surface 
structure under the working assumption that the star has a bare quark surface with an electrosphere, 
since electromagnetic effects dominate over QCD effects. A clear and necessary improvement is the 
introduction of surface magnetic fields. Neutron stars are born with intense magnetic fields of order 10^^ 
Gauss or higher, which may decay by over 5 orders of magnitude in the star's lifetime. If a neutron star 
undergoes a phase conversion to a bare quark star, the surface magnetic fields are expected to remain 
intact in geometry, unless interior magnetic fields are expelled by a superconducting state that displays 
the Meissner effect, in which case rapid magnetic reconnections occur, diminishing the field at surface 
to almost zero [19]. In this work, we will superimpose a static magnetic field B on the electrosphere 
of a quark star, and calculate the profile for the electrosphere as a first step toward characterizing the 
spectrum of a magnetized quark star. We expect considerable differences from the B — Q case since 
magnetized atmospheres of neutron stars are already known to be very different from non-magnetized 
ones [20]. Further, the inclusion of magnetic fields is a step toward a more realistic treatment of the 
quark star spectrum. 

In the next section, we reintroduce the established picture of the electrosphere for the purpose of a 
self-contained presentation, and mention the modifications due to the strong magnetic field. In section 
III, we quantify the electrostatic potential and electron number density in the magnetized electrosphere, 
and state qualitative impacts on the electrospheric photon emission. We conclude in section IV with a 
general discussion on the formation of Quark stars, with particular attention to the attractive Quark-Nova 
scenario which may be connected to r-process nucleosynthesis. 



II. DEGENERATE ELECTRON GAS IN A STRONG MAGNETIC FIELD 

Recently, a few authors [21] pointed out that the deficit of (massive) strange quarks due to surface 
effects on the star can lead to a thin charged skin a few Fermis thick, which can lead to a jump in the 
electric field at surface. The local charge neutrality condition is 2/3n„ — l/Sn^ — l/iug — ng = 0, with 
the appropriate low temperature (T ^ /i^) expressions in the bulk: 

rii^ — ^ = — ""- TT^ ; ■me<^mu,d^ms < iiq'iq . (1) 

Imposing the conditions of /3-equilibrium and baryon number conservation shows that /ie — m^/4/i to 
better than 1%. At surface, s-quarks are depleted due to surface tension effects by an amount [22] 
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leading to a charged layer at the surface of a bare quark star characterized by a surface charge density 
a — — (l/3)e(5ns and a discontinuous normal electric field with AE = Aira. For "cold" electrospheres, 
corresponding to the temperature bound T < 5 x lO^^K, the electric field is directed inward just under- 
neath the surface and outward just above it. The electrostatic potential and the electron number density 



are substantially increased of. the case when surface effects are ignored (NB. The corresponding surface 
depletion of light quarks is negligible.) [21]. 

The large electric field leads to the Schwinger instability of the QED vacuum [23] , resulting in electron- 
positron pair-emission, which can be an additional source of photons in the electrosphere, which normally 
radiates photons via 2^- 3 QED processes. A model cooling calculation of bare strange stars with 
electrospheric effects and the inclusion of color superconductivity (which can alter the specific heat 
and thermal conductivity of quark matter) has also been performed [13]. The conclusion is that bare 
strange stars will display Supcr-Eddington photon luminosities at surface temperatures T > lO^K, with 
a spectrum that is hard enough to distinguish it from thermally radiating neutron stars. At lower 
temperatures 6 x 10^K< T < lO^K, the bulk of the luminosity comes from electron-positron pairs which 
subsequently imprint a wide annihilation line on the non-thermal spectrum. Thermal emission is much 
suppressed owing to plasma frequency effects of quark matter and the electrosphere [24] . Further cooling 
to temperatures T < lO^K results in a non-thermal spectrum dominated by bremsstrahlung photons from 
electron-electron collisions in the electrosphere [14]. 

As mentioned in the introduction, a magnetic field is expected to exist at surface, altering the elec- 
trosphere profile and its emission properties. The electron motion perpendicular to the magnetic field 
is quantized into Landau orbitals enumerated by the quantum number ul = 0, 1, 2.. . For a relativistic 
fermion of charge q and mass m, such as the quarks and electrons in our case, the energy in a z-directed 
magnetic field B is given by 

e„ = Jm'^c'^ + c^pI + 2mc'^huJBn = mc^E; 2n = 2{nL + - + a); ujb = — , (3) 

2 mc 

where a = ±1/2 is the spin-degeneracy of the fermion. The n — state is non-degenerate while n > 
states are doubly degenerate. We express the magnetic field strength in convenient units '&' such that 

b = TT = 75 — for electrons (4) 

and take B = 4.414 x lO^^G (corresponding to 6 = 1) as an upper limit on the surface magnetic field 
of quark stars. Although higher fields may exist in magnetars and AXPs (anomalous X-ray pulsars), 
their interior magnetic fields would likely be large enough to inhibit the formation of quark matter [25]. 
For a given magnetic field, both temperature and matter density determine the filling fraction of the 
Landau levels, and hence the degree of quantization. From Eqn.(3), the highest Landau level occupied 
by a fermion with energy e„ is 

. ^ ^' - 1 
flL,max : "max = Int(l/); 1/ = — . (5) 

Since the quarks are highly relativistic, E ^ 1 and many Landau levels are occupied for the light 
quarks, and even more for the strange quark, so that quantization effects may be ignored. Further, 
when the temperature T ^ Tg where Tb = hLOs/kB {kB is the Boltzmann constant), thermal effects 
smear out adjacent Landau levels, and the magnetic field is effectively non-quantizing, irrespective of 
the degeneracy. This effect is more important for quarks, which have a larger mass, and hence smaller 
cyclotron frequency, than electrons. We thus expect no modifications to the quark distribution due to 
magnetic fields of typical magnitudes considered in this work. 

For electrons, quantization effects become important as we move outward from the electrosphere, and 
more electrons reside in the lower Landau levels. The number density of electrons in a magnetic field is 
given by 



For a cold degenerate electron gas which satisfies T ^ fXg, we may approximate 
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From Eqn.(7), the density at which the second Landau level becomes populated is 

^ (jUe^B?!^ ^ ^24 X 103"63/2 cm-3 . 
This is already less than the typical electron density at surface in the B 
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case, where Uf. ^ 
10"'' — 10""^/^ [26]. Therefore, several Landau levels will be occupied, and the i?-field is non-quantizing. 
In fact, for i^ > 1 and T ^ ^j.^, it follows from Eqn.(7) that [27] 
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which is the usual expression relating the number density and chemical potential for a relativistic 
degenerate Fermi gas in the B = case. For the range of temperatures and densities considered in 
this work, T <C Tb = hLOs/kB holds, so that temperature effects on electron occupation numbers is 
negligible. The situation is summarized in Fig. 1, which shows how the number density of electrons, as 
obtained from Eqn.(6) behaves as a function of electron chemical potential, temperature and magnetic 
field. The truncation of the curves at z^ = for finite T represent the fact that for lower (negative) 
chemical potentials, and hence densities, all electrons reside in the n = state. We note that only for 
temperatures T > Tb does the number density deviate significantly from the T — curve (NB. This is 
due to the smearing of the Landau levels by temperature fluctuations.). We therefore continue to use the 
zero temperature approximation to the electron density in the rest of the paper. 
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FIG. 1. The number density of electrons (scaled to nl) from Eqn. (6) as a function of the filling factor u for 
various temperatures. The solid curves are for non-zero magnetic field, and display a loss of the Landau level 
structure with increasing temperature. 



III. ELECTROSPHERE PROFILE WITH A MAGNETIC FIELD 



We would now like to obtain the electron density profile as a function of distance in the electrosphere. 
This will be modified from previous results [11,24] which were in the absence of a magnetic field. The 
baryon density at the surface of a strange star, where the quark pressure vanishes, is expected to be in the 
range (2 — 5)po, corresponding to quark chemical potentials in the range 300-400 MeV. At such densities, 
the strange quark current mass is not ignorable, with rus ~ 150 MeV. Local charge neutrality then implies 
an electron chemical potential at surface ^e ~ 15 — 20 MeV. This picture of a homogeneous surface as 
opposed to a crusty one with quark nuggets is justified by the large value of the surface tension given 



our parameter values [28] . The surface depletion of strange quarks at the surface results in a positively 
charged layer, which is two orders of magnitude thinner than the electrosphere. For cold electrospheres, 
this charged layer results in oppositely directed electric fields at the star's surface as well as an enhanced 
electron number density rie 3> ?t.*. 

For a general magnetic field with poloidal and toroidal components at the surface, the electron distri- 
bution can, in principle, be obtained using the formalism of magnetic field density functional theory or 
by Green function techniques, both of which are complicated tasks. As a first study, we adopt a simple 
picture of a constant magnetic field directed normally to the surface and ignore Coulomb interactions and 
exchange effects, which are justified approximations given the extreme degeneracy being explored in the 
problem. We employ the self-consistent Thomas-Fermi approximation and present a numerical solution 
of the Poisson equation 

V20(f,B) = -p(f,B)=ene(f,B). (10) 

In the plane parallel approximation B = Bz, f — z, and equilibrium requires /ie(fie('Z)) ~ e0(z) = 0. 
For arbitrary magnetic fields, the equilibrium constraint involves magnetic forces -Fmag = J x -B 7^ where 
j is the electron current density. With our approximations, the problem is an electrostatic one, and the 
electrostatic potential outside (z > 0) the star solves 
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In the limit iV 3> 1, when the electrons populate several levels, we recover the expression for the 
non-magnetic case 

V20+ = ^eV+ , (13) 

while for A^ = 0, we obtain 

V^0, = ^^^ , (14) 

the expression which has been derived previously for the non-relativistic case [29] but applies equally 
well to the relativistic case as long as only the 0**^ level is populated. It is interesting to note that the 
above equation leads to an exponentially decaying electron profile at large distances from the surface, 
contrary to the 1/z fall-off in the non-magnetic case. As r^-i- — > 1, we expect interactions and exchange 
effects to be important, but since this is only a minute region of the electrosphere, we do not correct for 
it. 

Eq. (11) can be solved numerically once the boundary conditions are determined. Clearly, ry_(_ ^ 1 as 
z ^ 00 (this corresponds to vanishing electron number density and is numerically implemented by the 
shooting method to obtain the correct asymptote at sufficiently large z. This method is also practically 
applicable to the case of suspended crusts above a quark star; see Ref. [16].). The second boundary 
condition is determined at the surface z = 0, where (f> = (f>o ^ t] = tjo and ijo must be determined by the 
discontinuity in the electric field V(/)+ — V(/)_ = Ana with (j>- = (j>{z < 0). The electrostatic potential </>- 
inside the star is almost identical to the non-magnetic case since the electron density deep inside the star, 
given by fie = ni^/Ang ~ 10 MeV ^ irieC^ — 0.5 MeV so that for b < 10.0 fields, quantization effects are 
negligible. In this case, (/>_ — > 0^ as z — > —00 where the relation between 0o and (j)q is given by 

4^'^ = '5(£) ['^o + (09-0o){('^o + ^g)' + 202|]i/2 ^ ^^5) 

wherein the prefactor corrects an error in equation (30) of Ref. [21]. For our choice of ile = 11.25 
MeV (corresponding to m^ = 150 MeV and ^, = 500 MeV), </ig = 37.4 MeV, 0o = 578 MeV and hence 
?7o = 341 MeV. The boundary conditions having been specified, the results for the electrostatic potential 
(j)+, (j)^ as obtained from Eq.(lO) are illustrated in Fig. 2, and compared to the corresponding expressions 
for the non-magnetized electrosphere with the same boundary conditions. 
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FIG. 2. Left panel: The electrostatic energy (scaled to rrisC^) at the edge of the star, from Eqii. (10) for the 
unmagnetized {b = 0) and magnetized (6 = 1.0) cases. Right panel: The corresponding electron number density 
(scaled to n*) at the surface. 

We observe that while the electrostatic potential inside the star is unaffected by magnetic field strengths 
of 6 '^ 0{1), that outside the star is increased in comparison to the unmagnetized case. A similar 
finding for the number density reflects the fact that the electrons become localized in configuration space, 
their wave functions compressed into the Larmor radius, so that the Pauli repulsion is decreased and 
electrons can be squeezed together in greater numbers than in the unmagnetized case. This effect is most 
pronounced when the quantization effects are large, as can be seen by the wiggles that demarcate the 
filling of the Landau levels (right panel of Fig. 2). In both figures above, the curves meet at the surface 
(z = 0) where they satisfy the same boundary condition. At distances far from the surface, the curves 
cross each other since the total number density is the same in both cases. The number density falls off 
exponentially for the magnetized case while it falls off as 1/z for the unmagnetized case. 

In the absence of a thin (lO^^Mo) crustal layer of accreted matter above the strange star, the spec- 
tral properties of a strange star are linked to radiative transport in the electrosphere. This is because 
radiation from the stellar interior is cutoff by the high plasma frequency (wp ~ 20 MeV) of dense quark 
matter, and Landau-Pomeranchuk-Migdal (LPM) suppression of bremsstrahlung photons from quark- 
quark collisions [24]. For an unmagnetized electrosphere, the emission at high luminosities is dominated 
by e+e~ pair emission so long as the temperature is high enough. At lower luminosities, the emission is 
dominated by photon bremsstrahlung from electron-electron collisions in the electrosphere. The spectrum 
is non-thermal due to the finite plasma frequency of the degenerate electron gas and large mean free path 
for photons in the electrosphere. 

The effect of the intense magnetic field on electrospheric emission is currently being investigated and 
will be reported elsewhere; here, we only state some qualitative effects. The plasma frequency plays 
a vital role in the emission properties of the electrosphere, since lower frequency photons are rapidly 
damped by the complexion of the refractive index 
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This collective effect of the plasma is important at high electron density (inner regions of the electro- 
sphere) where the Debye length is smaller than the mean free path. Note that despite the large electron 
density, the large degeneracy limits the available scattering states of electrons, leading to a large mean 
free path for photons. In the presence of a magnetic field, when quantization effects become important, 
collective effects are determined by the Larmor radius rather than the Debye length, and the reflectivity 
for normal incidence becomes [30] 
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When Wp <C w <C lob , the reflectivity is considerably diminished from the unmagnetized case, and the 
electron plasma will radiate more strongly. Further, the dielectric tensor will modify the Poisson equation, 
resulting in the breakdown of the plane parallel approximation employed here. Thus, our estimate of 
the profile is simplistic and must be made self-consistent. The intense magnetic field will also lead to 
polarization effects. Taking into account the thermal conductivity of the electrons in a strong magnetic 
field, a non-uniform surface temperature distribution is expected. The inclusion of toroidal components 
of the magnetic field can smoothcn out some of these non-uniformities, and has been linked to the largely 
uniform blackbody spectrum of some isolated neutron stars [31]. It would be interesting to see if the 
electrospheric surface of the strange star displays similar spectral features. In particular, the possibility 
of having bound positronium in the magnetized electrosphere should be investigated. 

IV. STRANGE STARS AND ASTROPHYSICS: THE QUARK-NOVA AND THE i?-PROCESS 

Strange quark stars have been linked to some of the most energetic and violent phenomena in astro- 
physics, such as Gamma-ray bursts (GRBs) [32] , type I X-ray superbursts [33] and Soft-gamma repeaters 
(SGRs) [34] in what is being perceived as a new paradigm in astrophysics - the connection of high-energy 
astrophysics with exotic forms of dense matter. The energy released in the phase conversion of a hadronic 
star to a quark star is typically about 10^^ ergs, similar in magnitude to the energy of observed gamma- 
ray bursts, and accretion and flash heating on a strange star surface is a viable model for SGR flares 
and superbursts. Further, the emission spectra of the long bursts can also be explained by invoking the 
presence of strange quark matter and its subsequent transition to a color-superconducting state [35] . All 
of these scenarios are predicated upon the existence of absolutely stable three-flavor quark matter at high 
density. 

Here, we discuss the formation of a quark star in a Quark-Nova, a scenario in which nuclear matter 
at high density in the core of a neutron star undergoes a phase transition to two flavor (up and down) 
quark matter, which then undergoes equilibrating weak reactions to form three flavor quark matter. 
The advantage is that this docs not require the improbable simultaneous and spontaneous conversion 
to strange matter (via weak interactions) of several neutrons within a small volume. In addition, the 
Coulomb barrier-free absorption of neutrons can enlarge the quark phase, so that the whole neutron star 
is essentially converted to strange quark matter, thus forming a quark star. The dynamics of this process 
is only beginning to be explored [36] , and it is likely that some mass ejection from the surface of a neutron 
star can take place, aided by neutrinos [37] but powered mostly by pre-compression shock waves from 
the rapidly advancing quark conversion front [38] . 

There are various ways in which this conversion can be triggered. A rapidly rotating neutron star 
is slowed down principally by magnetic braking (and additionally by energy loss through gravitational 
waves), thereby reducing the centrifugal forces and increasing central pressures. The timescales and 
dynamics of this spin-down have been discussed in [39]. The probability of triggering a first order 
transition by forming a tiny lump of two flavor quark matter (quantum nuclcation) is exponentially 
sensitive to the central pressure [40]. This makes the transition much more likely as the star spins down. 
Material accumulated from a fallback disk around a cooling neutron star or accretion from a companion 
star, for example, in a low-mass X-ray binary system can also force the conversion. Even the Bondi 
accretion rate of 10^^ g/sec from the interstellar medium can lead to a mass increase of up to O.IMq 
within a million years which can be sufficient to compress the star beyond the minimum density for 
the phase transition. Other possibilities include the clustering of A-baryons at high density and seeding 
by energetic cosmic neutrinos or strangelets [11]. Deconfinement of quarks can also occur as early as 
the protoneutron star stage following a core-collapse supernova explosion if the central density is large 
enough. However, in this case, neutrino trapping in the hot and dense interior can push the transition 
density higher, delaying the collapse to the point where a black hole might be formed instead [41]. 

The amount of mass ejected and the dynamics of the Quark-nova make it a potential candidate for 
r-process nucleosynthesis. The decompressing neutron matter at the surface of the neutron star which is 
undergoing the phase conversion to a quark star will expand rapidly and chunks of neutron matter will 
fall apart, allowing beta-decays to heat the material to close to r-process like temperatures of T ^ 0.1 
MeV. Following the expansion of this dense cjccta (including /3-decays) down to neutron drip density 
using network calculations of Meyer [42] , we find that the matter ends up distributed over a wide range 
of elements from Z ^ 40 — 70. At this stage, the full r-process reaction network can be coupled to 



obtain the final abundance distribution, as well as the evolution of temperature, entropy and electron 
fraction. The most important parameter that governs the evolution of the above physical quantities is 
the expansion timescale of the ejecta, which is the time taken to evolve from a density of p ~ lO-'^^g/cc 
to p ~ lO^-'^g/cc. Assuming that the ejecta is on an escape trajectory and expands as a multiple of the 
free-fall timescale in a spherically symmetric, non-interacting manner, we obtain a timescale on the order 
of 0.01-1 milliseconds. Fig. 3 shows the r-abundance of the elements (labeled "network") from a Quark- 
Nova compared to the normalized solar r-abundance pattern. The overabundance for elements with mass 
number A > 230 is due to the omission of alpha-decays in the network. It is also notable that not all 
neutrons are consumed in this scenario, and that several light elements are produced with significant 
abundance. It is likely that Quark-Novae may be an additional contributor to the r-process material in 
the universe, presenting another link between dense quark matter and astrophysics. A detailed analysis 
of r-process nucleosynthesis from Quark-Novae, including contrasts to standard or primary astrophysical 
sites such as Type II supernovae and neutron star mergers, and comparisons to trends of heavy-element 
abundance ratios in metal-poor stars will be presented in a forthcoming article [38]. 
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FIG. 3. Abundance pattern of r-process elements from a Quark-Nova: The initial electron fraction Ye = 0.03, 
entropy per nucleon S/kn = 1 , temperature Tg — 1.0 and expansion timescale r = 1.0 millisecond. The scaled 
comparison to solar abundance is also shown. 



While several observable characteristics of strange stars have been conjectured in the literature [43,44], 
such as period-clustering, anomalies in braking index, global properties (mass and radius), and small 
spin-periods, spectral differences between quark stars and neutron stars have only been touched upon. 
If strange stars are made in Quark-Novae, and are linked to r-processing of the heavy elements, one 
can expect a most direct evidence of a quark star to come from 7-decays of unstable nuclei produced in 
the r-process. We expect the Quark-Nova ejecta to achieve 7-ray transparency sooner than supernova 
ejecta since Quark star progenitors (i.e., neutron stars) lack extended atmospheres, so that r-process only 
gamma emitters with lifetimes of the order of years (such as ^^^Cs, ^'^'^Ce, ^^^Eu and ^^^Os) can be used 
as tags for the Quark-Nova. 

In conclusion, this is an exciting time for quark star research. Recent developments in the theory of 
high density quark matter have provided novel ways of explaining the most energetic astrophysical events, 
and spurred re-examinations of old paradigms of strange stars. With increasing data now coming from 
myriad sources (egs. neutron stars that push theoretical limits on spin periods and radius; GRB spectra 
that link to quark stars as inner engines; abundance trends of heavy elements in ultra metal-poor stars 
that point to multiple sources of the r-process), it appears that quarks may play an important role in 
uncovering the secrets of many astrophysical puzzles. 
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